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SUMMARY  AND  CONCLUSIONS 

Waveform  optimization  is  a  technique  that  can  be  used  to 
adjust  any  or  all  of  the  radar  transmitter  degrees  of  freedom 
(DoFs)  based  on  observations  of  a  dynamic  environment 
caused  by  intentional  and/or  unintentional  interference  sources. 
These  DoFs  include  but  are  not  limited  to  pulse  shape,  pulse 
repetition  frequency  (PRF),  bandwidth,  and  operating  fre¬ 
quency.  Waveform  optimization  techniques  can  potentially 
provide  an  additional  tool  to  improve  radar  system  perfor¬ 
mance  in  severe  operating  environments.  This  paper  addresses 
the  use  of  waveform  optimization  in  radio  frequency  interfer¬ 
ence  (RFI)  environments,  particularly  in  the  situation  where 
the  RFI  signal  is  partial  band  or  is  scattered  into  the  radar 
receiver  main  beam  by  terrain  (hot  clutter). 

1.0  INTRODUCTION 

Waveform  optimization  is  a  technique  that  can  be  used  to 
adjust  any  or  all  of  the  radar  transmitter  degrees  of  freedom 
(DoFs)  based  on  observations  of  a  dynamic  environment 
caused  by  intentional  and/or  unintentional  interference  sources. 
These  DoFs  include  but  are  not  limited  to  pulse  shape,  pulse 
repetition  frequency  (PRF),  bandwidth,  and  operating  fre¬ 
quency.  While  adaptivity  on  receive  has  been  exploited  for 
quite  some  time  [1-3]  waveform  optimization  provides  the 
promise  of  the  development  of  a  fully  adaptive  radar.  Recent 
work  has  resulted  in  various  techniques  for  applying  this  opti¬ 
mization  [4-6].  Waveform  optimization  techniques  can  poten¬ 
tially  provide  an  additional  tool  to  improve  system 
performance  in  severe  operating  environments. 

This  paper  addresses  the  use  of  waveform  optimization  in 
radio  frequency  interference  (RFI)  environments,  particularly 
in  the  situation  where  the  RFI  signal  is  partial  band  or  is  scat¬ 
tered  into  the  radar  receiver  main  beam  by  terrain  (hot  clutter). 
There  can  be  significant  coverage  losses  due  to  terrain-scat¬ 
tered  interference.  For  finite  aperture  the  only  solution  may  be 
to  “bum  through”  which  will  impact  the  radar  time  line.  Wave¬ 
form  optimization  can  potentially  provide  an  alternative  to  this 
by  emphasizing  (matching)  to  the  target  and  de-emphasizing 
(anti-matching)  to  the  interference. 

Section  2  summarizes  the  application  of  waveform  optimi¬ 
zation  theory  to  RFI  mitigation  and  develops  a  bound  on  the 
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achievable  gain  of  the  optimal  waveform  in  terms  of  the  eigen¬ 
values  of  the  covariance  matrix  of  the  REI  signal.  Section  3 
summarizes  the  waveform  optimization  algorithms  employed 
in  the  study.  Section  4  provides  analysis  of  the  performance  of 
waveform  optimization  techniques  in  the  presence  of  band- 
limited  white-noise  REI  and  terrain-scattered  interference. 
Summary  and  conclusions  are  provided  in  Section  5. 

2.0  SUMMARY  OF  THEORY  AND  PERFORMANCE  BOUND 
We  begin  by  considering  the  optimum  receiver  configura¬ 
tion  for  a  signal  corrupted  by  wide-sense  stationary  (w.s.s) 
additive  colored  noise  (ACN)  and  follow  the  development  of 
[4-6].  This  configuration  is  shown  in  Eigure  1.  It  is  well  known 
that  the  signal-to-interference-plus-noise  ratio  (SINR)  is  maxi¬ 
mized  when  the  received  signal  x{t)  is  first  whitened  to  cancel 
the  ACN  and  then  detected  using  a  filter  matched  to  the  ‘whit¬ 
ened’  waveform  [8]  i.e., 

=  s(0*M0;  Vf(0  =  [Sv,(i-D]*. 
where  *  denotes  linear  convolution  and  superscript  *  denotes 
the  complex  conjugate.  Thus  the  output  is  then 

>’(0  = 

where  /z^(0  =  h^{t)  *  ■ 


n-jit) 

additive 
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Eig.  1.  Optimal  receiver  configuration  for  detecting  a  sig¬ 
nal  in  additive  colored  noise. 

Given  this  receiver  structure,  our  task  is  now  to  exploit  the 
degrees  of  freedom  (DoEs)  available  in  the  transmit  waveform 
to  maximize  the  SINR.  Again,  following  [4-6],  we  first  con¬ 
sider  the  SINR  of  the  optimally  filtered  received  signal: 

U 

SINRo  =  \\\s^{tfdt 

^  1 

2 

where  is  the  power  spectral  density  of  the  whitened  ACN 
and  Tj  and  Tj  are  the  limits  of  the  observation  interval.  To 
maximize  the  SINR,  we  must  maximize  the  energy  in  the  whit¬ 
ened  signal  over  all  possible  transmit  waveforms.  Thus  we 
wish  to  solve: 
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max  [  fs^(Ti)/z^(f-Ti)£/Ti 
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r,  VO  V  Vo  / 

and  we  have  made  the  assumption  that  the  transmit  waveform 
is  of  finite  duration  over  the  interval  of  [0,  7^  .  Rearranging 
terms  and  changing  the  order  of  integration  yields: 

T  T 


max  fs^(Ti)fiy(T2)/i^(f-Ti)  : 

Sj^ty 


0  0 

/r. 


jK(t-x^)h^(t-X2)dt 


Vr, 
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where  the  last  term  is  the  conjugate  of  the  kernel  (which  bears 
a  resemblance  to  the  autocorrelation)  function  of  h^{t)  that  we 
denote  K{x^,  T2) .  If  we  apply  Schwarz’s  inequality,  we  find 
the  optimal  waveform  SQ(t)  to  be  the  solution  of 

T 


that  is,  the  optimal  waveform  is  an  eigenfunction  associated 
with  the  maximum  eigenvalue  of  the  kernel  /f(Tj,  T2)  • 

If  we  consider  a  practical  implementation  using  discrete¬ 
time  signals,  we  have  the  receiver  structure  shown  in  Figure  2. 
Now  the  optimum  SINR  is  given  by 

II 2 


SINR 


2  //  H 
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where  H  denotes  the  Hermitian  transpose,  H*!!  is  the  L-2  norm, 
and  again  this  quantity  can  by  maximized  by  appropriate  selec¬ 
tion  of  Sy.  This  will  be  the  eigenvector  associated  with  the 

^  H 

maximum  eigenvalue  of  .  It  is  well-known  that  the 

whitening  filter  to  cancel  the  ACN  n  is  where 

=  ^[nn  ]  and  £[•]  is  the  expectation  operator.  Thus  R^, 
is  a  Toeplitz  matrix  with  the  appropriate  autocovariance  lags 
along  each  diagonal.  (This  can  be  verified  by  considering 
£’[({Rj,^^^}^n)({Rj,'  ^}^n)^]).  Thus  the  optimal  transmit 
waveform  is  Sq  =  where  is  the  eigenvector  associ¬ 


~min  ' 


ated  with  the  minimum  eigenvalue 

We  note  that  the  output  of  the  receiver  of  Figure  1  can  be 
written  as 


^TT 


St  =  Umin  ■  Since 


y  - 

When  the  ogtimal  waveform  is  used, 

R„  =  (H„H")  ^  can  be  written  as  UAU^,  where  U  is  a 

n  v  vy  Vv  2  ’ 

matrix  whose  columns  are  the  eigenvectors  of  Rj,  and  A  is  a 
matrix  whose  diagonal  elements  are  the  eigenvalues  of  Rj, ,  we 
can  write 


y 


UmmUA  U  X 


.  -1  H 


n 
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Fig.  2.  Discrete-time  form  optimal  receiver  configuration 
for  detecting  a  signal  in  additive  colored  noise.  Vectors  of 
the  various  time  series  are  indicated  by  Sj,  n,  x,  and 
^MF  while  represents  a  matrix  transformation. 


Thus  when  the  optimal  waveform  is  used,  the  whitening  filter 
stage  of  Figure  1  can  be  eliminated. 

If  we  define  and  as  the  minimum  and  maxi¬ 

mum  eigenvalues  of  R^j ,  respectively,  then  the  minimum  and 
maximum  eigenvalues  of  R^  are 

respectively.  We  then  see  that  the  SINR  is  bounded  as 
shown  [7]: 


^<SINRo<^. 

^max  ^min 

This  double  sided  inequality  also  leads  to  an  upper  bound 
on  the  achievable  gain  of  the  optimal  waveform  over  any  other 
arbitrary  waveform.  This  is  [7]: 

2 


K 


(1) 


where  g  is  the  relative  gain  between  a  radar  that  uses  the  opti¬ 
mal  waveform  and  one  that  uses  some  other  suboptimal  wave¬ 
form.  Therefore,  the  achievable  gain  is  characterized  by  the 
eigenvalue  spread  of  the  interference  covariance  matrix.  It  is 
interesting  to  note  that  when  the  radar  observation  interval  is 
long  the  eigenvectors  (eigenfunctions  in  the  continuous-time 
case)  become  complex  exponentials  (see  [8]  p.  205)  and  the 
bound  given  in  (1)  can  be  expressed  as  the  ratio  of  the  maxi¬ 
mum  and  minimum  value  of  the  interference  power  spectral 
density  (PSD)  over  the  system  bandwidth.  That  is. 


max  A(/) 
^“minVC/) 

where  N(f)  is  the  PSD  of  the  interference  signal. 


3.0  ALGORITHMS 

We  consider  three  methods  for  calculating  the  transmit 
waveform  based  on  the  interference  (temporal)  covariance. 
‘Method  1’  is  to  simply  chose  the  eigenvector  associated  with 
the  minimum  eigenvalue  of  the  interference  covariance  [4-6]. 
The  rational  is  to  place  the  transmit  energy  into  the  dimension 
with  the  least  amount  of  interference.  In  this  case,  since  eigen¬ 
vectors  have  a  unity  norm,  the  optimum  waveform  Sq  = 
is  normalized  to  unit  energy. 

‘Method  2’  takes  advantage  of  the  rank-deficient  nature 
the  interference  of  Figure  3  and  uses  the  sum  all  of  the  eigen¬ 
vectors  corresponding  to  eigenvectors  at  or  below  the  noise 
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level,  and  then  normalizing  the  result  to  unit  energy.  Thus  we 
have; 

K' 

So  =  Y  X  “-t 

k  =  1 

where  VLj^,k  =  are  the  ‘noise’  eigenvectors  (i.e.,  the 

eigenvectors  not  associated  with  the  ‘colored’  interference) 
and  y  is  chosen  so  the  result  has  unit  energy.  This  method  adds 
flexibility  for  the  optimum  waveform  to  spread  out  over  the 
entire  fraction  of  the  band  not  occupied  by  the  interference. 

‘Method  3’,  which  is  similar  to  the  approach  taken  in  [9], 
recognizes  that  traditional  waveforms  have  many  desirable 
characteristics  (constant  modulus,  ability  to  control  range  side- 
lobes  in  radar  applications,  etc.).  Thus  the  objective  is  to  find  a 
set  of  weights  such  that 

K' 

So  =  Y  X 

k  =  1 

is  ‘closest’  in  a  least-squares  sense  to  some  desired  waveform. 
The  constant  y  again  is  chosen  to  normalize  the  waveform  to 
unit  energy.  To  obtain  the  desired  result,  we  simply  need  to 
solve  the  equation: 

K' 

y  =  Z  =  UjT'a  =  Srf  (3) 

k  =  1 

where  =  [uj  Uj  ...  ;  tx  =  a2  ...  aj  ;  and 

is  a  vector  representation  of  the  desired  waveform  (below 
we  provide  examples  using  a  linear  frequency  modulated 
(LFM)  waveform).  To  solve  for  a  ,  we  simply  multiply  the  last 
two  expressions  of  (3)  by  U^,  and  exploit  the  unitary  nature  of 
.  Thus  we  find  the  result 

So  =  yU^,u"s^.  (4) 

4.0  SIMULATIONS 

Each  of  the  techniques  discussed  in  Section  3  was  ana¬ 
lyzed  in  terms  of  SINR  performance  and  ability  to  maintain 
desirable  radar  waveform  characteristics  such  as  constant  mod¬ 
ulus  and  low  range  sidelobes. 

4. 1 .  Band-Limited  Interference 

In  the  presence  of  interference  that  is  white  noise  that  has 
been  band-limited,  the  performance  gain  will  come  principally 
from  isolating  the  transmit  energy  into  interference-free 
regions  of  the  system  band.  Thus  these  gains  will  of  course  be 
greatest  when  the  interference  bandwidth  is  less  than  the  sys¬ 
tem  bandwidth.  To  examine  this,  we  consider  the  scenario 
illustrated  in  the  Figure  3,  where  the  interference  is  shown  to 
be  some  fraction  of  the  system  bandwidth.  While  we  do  not 
expect  a  difference  in  performance  based  on  the  location  of  the 
interference  within  the  system  band,  we  consider  both  the  off¬ 
set  and  centered  cases  shown  to  verify  this  assumption. 

We  first  consider  the  case  of  a  1  MHz  system  band¬ 
width  and  a  100  kHz  interference  bandwidth.  In  Figure  4, 


Fig.  3.  Interference  spectra  vs.  system  bandwidth. 

the  spectra  of  the  interference  and  the  optimal  waveform 
(Method  1)  are  shown  for  an  offset  case,  along  with  the  pulse- 
compressed  (detected)  waveform  moduli  as  a  function  of  time. 
The  results  for  Method  2  for  the  same  case  are  shown  in  Figure 
5.  Here  we  see  a  result  similar  to  Method  1. 


Power  Spectra 


Pulse-Compressed  Time  Series 


time  (ms) 


Fig.  4.  Top:  waveform  spectra.  Method  1  optimization. 

Bottom:  corresponding  time  series. 

This  is  in  contrast  with  Method  3,  for  which  results  are 
shown  in  Figure  6.  Here  we  see  that  the  pulse-compressed  out¬ 
put  is  much  more  well  behaved  and  retains  much  of  the  range 
sidelobe  characteristic  of  the  LFM  waveform. 

A  summary  of  SINR  performance  is  shown  in  Figure  7. 
We  consider  gain  in  SINR  relative  to  a  LFM  waveform.  We  see 
that  none  of  the  three  methods  attains  the  bound  when  the 
interference  is  a  smaller  fraction  of  the  system  bandwidth, 
although  substantial  gains  do  result.  In  these  particular  cases, 
large  gains  are  possible  as  a  result  of  concentrating  all  of  the 
interference  energy  in  the  ‘clear’  part  of  the  spectrum.  (In  other 
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Power  Spectra 


Pulse-Compressed  Time  Series 


time  (ms) 


Fig.  5.  Top:  waveform  spectra,  Method  2  optimization. 
Bottom:  corresponding  time  series. 


Power  Spectra 


Pulse-Compressed  Time  Series 


Fig.  6.  Top:  waveform  spectra,  Method  3  optimization. 
Bottom:  corresponding  time  series. 


scenarios,  this  may  not  be  possible,  and  smaller  gains  may 
result.)  Thus,  considering  the  desirable  sidelobe  and  modulus 
characteristics  of  Method  3,  this  would  typically  be  the  pre¬ 
ferred  method  of  the  three  described  herein.  Finally,  we  also 
note  as  expected  that  there  is  little  SINR  performance  differ¬ 
ence  among  the  methods  or  between  the  ‘centered’  and  ‘offset 
cases. 


Performance  Gain,  Centered  Case 


Fig.  7.  SINR  performance  gain  summary.  Top:  ‘centered’ 
case.  Bottom:  ‘offset’  case.  Here,  ‘gain’  denotes  the  poten- 
tiaFachieved  gain  in  SINR  relative  to  a  LFM  waveform. 

We  note  that  the  performance  gains  shown  herein  are  rela¬ 
tive  to  the  use  of  an  LFM  waveform  without  the  use  of  the 
whitening  filter  of  Figure  1 .  Future  analyses  will  consider  the 
impact  of  the  whitening  filter  in  combination  with  optimal 
waveforms. 

4.2.  Terrain-Scattered  Interference 

Terrain-scattered  interference  (TSI)  [10-13]  results  when 
terrain-scattered  signals  from  an  interference  source  are  scat¬ 
tered  into  the  radar  receiver  main  beam.  Space-time  adaptive 
signal  processing  can  mitigate  the  interference  but  at  a  perfor¬ 
mance  cost  when  radar  clutter  and  TSI  must  be  jointly  miti¬ 
gated. 

Since  the  multiple  delay  (and  Doppler-shifted)  signals  of 
the  interference  that  comprise  the  TSI  will  be  colored,  the 
potential  exists  to  optimize  the  radar  transmit  waveform  to 
occupy  temporal  eigenspace  where  the  interference  is  weak, 
thereby  improving  the  radar  signal-to-interference  ratio  with¬ 
out  using  space-time  degrees  of  freedom.  Even  small  improve- 
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ments  can  be  significant  in  reducing  bum-through 
requirements  or  interference  rank.  However,  motion  of  the 
interference  and/or  receiver  platform  results  in  the  interference 
being  non-stationary  [10,  11],  thus  latency  becomes  an  issue. 

The  interference  signal  was  assumed  to  be  white  covering 
the  full  radar  bandwidth.  The  signal  as  received  on  a  direct 
path  is  not  colored,  and  if  not  spatially  notched,  diminishes  the 
advantage  of  the  optimized  waveform.  The  radar  adapts  to  the 
colored  interference  by  estimating  the  temporal  covariance 
matrix  from  the  received  TSI  and  then  transmitting  an  opti¬ 
mized  waveform.  The  TSI  covariance  is  not  stationary  as  a 
result  of  Doppler  associated  with  the  receiver  and  interference 
source;  consequently,  the  advantage  from  using  the  optimized 
waveform  decays  with  time. 

The  analysis  was  carried  out  using  ISL’s  SCATS  model 
[14-16]  to  characterize  the  TSI  in  terms  of  power,  delay,  and 
Doppler  spread  of  the  scattered  signals.  The  scenario  for  this 
analysis  includes  a  UHF  interference  source  at  7,400  m  altitude 
flying  due  north  at  a  speed  of  165  m/s  over  WSMR,  NM.  The 
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Fig.  8.  Terrain  scatter  power  calculated  using  SCATS. 

radar  is  flying  due  north  at  a  speed  of  125  m/s  207  km  to  the 
northeast  of  the  interference  source.  The  interference  (Tx)  and 
radar  (Rx)  locations  are  shown  in  Figure  8.  A  white  noise  inter¬ 
ference  signal  of  1  MHz  bandwidth  is  assumed  matched  to  the 
radar  bandwidth  of  1  MHz.  As  a  result,  this  represents  a  stress¬ 
ing  case  for  which  there  is  no  ‘empty’  portion  of  the  spectrum 
for  the  waveform  to  fill. 

SCATS  calculates  the  scatter  power  from  each  scattering 
cell  on  the  terrain.  For  this  analysis  a  random  phase  chosen 
uniformly  from  the  interval  [0,  27i]  was  assigned  to  each  cell 
since  the  path  length  differences  among  the  individual  scatter¬ 
ing  elements  that  make  up  the  cell  span  many  RF  wavelengths. 
The  scattered  signal  from  each  terrain  cell  is  then  a  complex 
number  whose  magnitude  is  the  square  root  of  the  power 
received  from  the  cell  and  whose  phase  is  the  assigned  random 
phase  modified  by  the  propagation  delay  of  the  bistatic  path. 
Methods  of  calculating  the  TSI  covariance  can  be  found  in 
[10,  11] 

Analysis  was  performed  by  optimizing  the  waveform 
based  on  the  interference  covariance  observed  on  a  relatively 
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Fig.  9.  Receive  antenna  pattern  at  0  degrees  elevation. 

narrow  spatial  receive  beam.  For  this  example,  the  antenna  is  a 
horizontal  array  of  50  dipoles  spaced  at  half  wavelength.  Bore- 
sight  of  the  array  is  at  270°  azimuth.  The  bearing  to  the  inter¬ 
ference  source  is  261°  azimuth  or  9°  from  boresight.  The  beam 
is  steered  to  265°,  which  is  just  4°  in  azimuth  from  the  interfer¬ 
ence.  All  angles  are  in  true  azimuth.  The  beam  pattern  is  shown 
in  Figure  9.  A  suppression  of  the  backlobe  by  20  dB  can  be 
observed  in  the  pattern. 

Gain  of  the  optimized  waveform  relative  to  an  LFM  is 
shown  in  Figure  10.  Here  the  direct  path  is  assumed  to  be  spa¬ 
tially  suppressed  and  the  coloring  of  the  spectrum  due  to  the 
interference  multipath  allows  for  some  gains  when  using  an 
optimized  waveform.  However,  the  non-stationarity  of  the 
interference  due  to  the  platform-motion  induced  Doppler  shifts 
results  in  degradation  with  time.  The  waveform  synthesized 
from  50  eigenvectors  has  decayed  to  4  dB  after  20  ms.  We 
would  expect  results  to  vary  with  scenario  and  varying  plat¬ 
form  speeds,  etc. 
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Fig.  10.  Gain  of  using  an  optimized  waveform  synthesized 
at  time  0  as  a  function  of  time  parametric  in  the  number  of 
eigenvalues  spanned  by  the  optimized  waveform  for  the 
steered  beam  (Method  3).  The  interference  covariance  was 
characterized  using  100  lags. 
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5.0  SUMMARY  AND  FUTURE  WORK 

The  radar  waveform  can  be  optimized  to  improve  the  sys¬ 
tem  SINR  performance  by  taking  advantage  of  spectrum  that  is 
not  occupied  by  the  interference.  Issues  involve  a  trade-off 
among  SINR  improvement,  radar  detection  (range  resolution 
and  sidelobes),  deviation  of  the  waveform  from  constant  mod¬ 
ulus,  additional  computational  requirements,  and  persistence  of 
the  improvement  over  time  until  the  waveform  is  recalculated. 
Even  a  few  dB  of  additional  interference  suppression  can  have 
a  substantial  impact  on  radar  performance  allowing  greater 
range,  reduced  power,  faster  surveillance  revisit  rates,  or 
reduced  ‘burn  through’  requirements  during  engagements. 

Simulations  show  that  in  the  presence  of  partial  band  inter¬ 
ference  an  optimized  waveform  can  achieve  substantial  gains 
(order  of  30  dB)  relative  to  a  conventional  LEM  waveform. 
The  gains  and  quality  of  the  resulting  waveform  are  a  function 
of  the  relative  bandwidths  of  the  interference  and  radar 
receiver. 

TSI  simulations  showed  that  terrain  scatter  colors  the 
interference  signal  spectrum.  Eor  the  site-specific  scenario  sim¬ 
ulated,  gains  in  excess  of  20  dB  signal-to-interference  ratio  can 
be  calculated  but  at  the  expense  of  desirable  properties  of  the 
radar  waveform.  When  the  waveform  is  constrained  to  behave 
more  like  a  conventional  LEM  waveform,  the  improvements 
are  on  the  order  of  5-6  dB.  Platform-motion-induced  nonsta- 
tionarity  of  the  interference  results  in  a  reduction  in  perfor¬ 
mance  with  time  relative  to  the  covariance  estimate. 

It  was  generally  found  that  constraints  are  necessary  on  the 
optimized  waveform  to  provide  usable  radar  performance. 
Eorcing  the  optimized  waveform  to  look  more  like  a  radar 
waveform  with  desirable  characteristics,  e.g.  LEM,  improved 
the  compressed  pulse  resolution  and  sidelobes  greatly.  There  is 
a  trade-off  between  the  signal  to  interference  improvement  and 
the  ‘quality’  of  the  radar  waveform.  As  the  eigenvector  sub¬ 
space  used  for  the  optimized  waveform  shrinks,  the  optimized 
waveform  exhibits  poorer  pulse  compression  and  unsteady 
amplitude 

Euture  work  will  consider  issues  regarding  the  impact  of 
the  use  of  a  whitening  filter  on  system  performance,  a  broader 
range  of  interference  scenarios,  covariance  estimation  issues, 
and  other  methods  for  constraining  the  waveform  to  have 
desirable  characteristics. 
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